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Abstract 

I summarize recent work on non-Fermi liquids within certain generalized 
Anderson impurity model as well as in the large dimensionality (D) limit 
of the two-band extended Hubbard model. The competition between local 
charge and spin fluctuations leads either to a Fermi liquid with renormalized 
quasiparticle excitations, or to non-Fermi liquids with spin-charge separation. 
These results provide new insights into the phenomenological similarities and 
differences between different correlated metals. While presenting these results, 
I outline a general strategy of local approach to non- Fermi liquids in correlated 
electron systems. 
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I. INTRODUCTION 



Since its introduction four decades ago, Landau's Fermi Liquid theory has been the 
standard model for interacting many-fermion systems0. The theory postulates that at low 
energies only the quasiparticle excitations play an essential role. The quasiparticles, essen- 
tially dressed fermions, can be adiabatically connected to certain non-interacting fermions 
as we turn off the interaction strength. The Fermi liquid description has been successful not 
only for the weakly interacting electrons in simple metals, but also for strongly correlated 
fermion systems. In this latter category are liquid ^Heii for which the Fermi liquid theory 
was initially formulated, and the metallic states of V0O3 based compound^!!, the prototype 
material displaying the Mott transition phenomenontl. Also included are the "conventional" 
heavy fermiongBi such as CeCug and UPts, in which the mass of the quasiparticles are 
enhanced by as much as hundreds from the band-theory predictions. 

In recent years, a number of strongly correlated materials have emerged which show 
physical properties anomalous in the context of the canonical Fermi liquid theory. These 
include, in addition to the much studied high Tc copper oxide superconductors!, a class of 
novel heavy fermions0, d— or f— electron based metals close to quantum criticalityS'lll, 
quasi-one-dimensional materialsl^llJ, as well as certain artificially fabricated metallic point 
contactsEl'0. 

The theoretical question, then, is: under what conditions do electron correlations lead to 
a breakdown of Fermi liquid theory? In the past few years, several theoretical approaches 
have been taken to address this question. One approach builds on our understanding of 
the breakdown of Fermi liquid theory in one dimension. For weakly interacting one dimen- 
sional fermion systems, the perturbative renormalization group (RG) leads to the g-ology 
classification of spatially homogeneous metallic states. The possible states are Luttinger 
liquids or those with divergent CDW, SDW or superconducting correlation functionsEJ'lia. 
In dimensions higher than one, perturbative RG analysis has shown that, the Fermi lifluid 
theory does describe weakly interacting fermion systems with a regular density of statestatl. 
The mechanism for the breakdown of Fermi liquid theory is necessarily non-perturbative in 
interaction strengthil. 

An alternative approach to non-Fermi liquids uses local physics as a starting point. The 
motivations behind this approach are multi-fold. First of all, most of the correlated elec- 
tron systems are transition-metal, rare earth or actinides based compounds. The dominant 
electron-electron interactions in these systems are local in space. This is the result of quan- 
tum chemistry: the partially-filled d— or /— orbitals are much more contracted than the 
s— and orbitals of the simple metals and covalent semiconductors, making the intra-site 
Coulomb interactions by far the largest interaction parameter. Secondly, Anderson- and 
Kondo- like impurity models have been studied extensively for more than three decadesil. 
In particular, the multi-channel Kondo problem has long been recognized to display RG 
fixed points of the non-Fermi liquid varietyS. Finally, the large D dynamical mean field 
theoryi'i opens the door for systematic treatment of the competition between local dy- 
namics and spatial fluctuations. 

The work summarized here covers a specific source of local physics towards non-Fermi 
liquids, namely the competition between local charge (valence) fluctuations and spin fluctu- 
ations. This belongs to the domain of mixed valence physics, a classic problem in condensed 
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matter theory. It is different from the multi-channel physics for non-Fermi liquiddHil'El. In 
the remainder of this section, I introduce the problem, define the models, and summarize the 
essential new results. More detailed discussions of the underlying physics are given in the 
subsequent sections: Sections II and III focus on the single impurity generalized Anderson 
model, and Sections IV and V discuss the two-band extended Hubbard lattice model. 



A. Phenomenological considerations 

Extensive studies on the heavy fermion metals have led to a canonical picture for the 
formation of a Fermi liquid in metals with strong local electron-electron interactions. Fig. 
1 illustrates the point. Plotted here are the temperature dependence of the Cu-site NMR 
relaxation ratec3 (1/Ti) and that of the electrical resistivityc2l in CeCue, one of the so-called 
"vegetable" heavy fermions. At asymptotically low temperatures, the NMR relaxation rate0 
is linear in temperature, while the electrical resistivity is quadratic in temperature. Both 
are characteristic of quasiparticle contributions. Simplistically speaking, the number of 
thermally excited spin excitations is proportional to temperature as a result of the Fermi- 
Dirac distribution. And each spin excitation contributes a temperature-independent term 
to the fiipping rate of nuclear spins, but a T— linear term to the quasiparticle scattering 
rate. The latter is again due to the Fermi statistics, which reduces the phase space for 
quasiparticle-quasiparticle scatterings. The experimental data behave very differently at 
temperatures above about 5 — lOK. Here, the NMR relaxation rate becomes essentially 
temperature-independent, while the temperature dependence of the electrical resistivity is 
characteristic of Ko„do-=catteri„g from lical momentfl. A local moment picture serves as 
a better starting point to describe the /—electron degrees of freedom in this temperature 
range. 

This crossover phenomenon provided the phenomenological basis for the canonical the- 
oretical picture for heavy fermion metals. In this picture, the /—electrons cross over from 
incoherent moments at high temperatures to being part of the renormalized quasiparticles 
in a coherent Fermi liquid at the lowest temperatures. This is the lattice analog of the broad 
crossover that is known in the solution to the single impurity Kondo problemE^E^I. In the 
single impurity problem, the characteristic crossover temperature is the Kondo energy. In 
the lattice case, the crossover temperature relates to the coherence energy scale below which 
the elementary excitations are quasiparticles with heavy mass. The coherence energy acts 
as the renormalized Fermi energy for the low energy quasiparticle excitations. 

This canonical picture appears to break down for a set of novel /—electron materials^. 
At low temperatures, these materials typically have an electrical resistivity linear in tem- 
perature, accompanied by anomalous features in a host of other physical properties. The 
precise mechanisms for these low temperature non-Fermi liquid phenomenologies are at this 
stage unclear. We refer the readers to the contribution of Miranda et al.0 in this volume 
for a survey of theoretical ideas. The crossover from high to low energies in these systems 
are only beginning to be addressed^!. 

In the normal state of the high Tc cuprates, the spin dynamics appear to show a crossover 
qualitatively similar to that of the heavv fermions. Shown in Fig. 2 are the temperature 
dependences of the NMR relaxation rateE^ and electrical resistivityCZl in the optimally doped 
La2_xSrxCu04. The temperature dependence of the NMR relaxation rate behaves in a way 
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reminiscent of that of CeCue- It has the asymptotic low temperature T— hnear behavior, 
crossing over to an essentially temperature-independent behavior at high temperatures. The 
crossover temperature scale is about 300K, much higher than that of the heavy fermions. 
The spin excitations can be thought of as the quasiparticle-quasihole continuum at low tem- 
peratures, but as excitations derived from local moments with short range antiferromagnetic 
coupling at high temperatures. This picture is corroborated by the neutron scattering re- 
sults. The low energy incommensurate peaks are most naturally accounted for in terms of a 
quasiparticle contribution, while the high energy broad background is naturally interpreted 
in terms of short range local moment correlationSH. 

However, when it comes to the temperature dependence of the electrical resistivity, the 
analogy with the heavy fermions stops. As can be seen in Fig. 2(b), the electrical resistivity 
is linear in temperature over essentially the entire temperature range. 

By now, there exist strong evidences that the dominant contribution to the electrical 
resistivity in the high cuprates comes from electron-electron scatteringSill, as is the case 
for the heavy fermions. The NMR relaxation rate is of course dominated by the electronic 
contributions. It is therefore quite unusual that electrons in these systems yield very similar 
magnetic responses (albeit with different energy scales), but entirely different charge trans- 
port properties. At the microscopic level, both the heavy fermions and copper oxides can 
be described by a model with a strongly correlated band and a weakly correlated one. The 
strongly correlated band is formed from the /— orbitals in the heavy fermions and from the 
3^2:2-2/2 orbitals in the cuprates. The weakly correlated band is from the non— / orbitals in 
the heavy fermions and from the oxygen 2p orbitals in the cuprates. For heavy fermions, the 
point of departure for most theoretical work is the Anderson lattice model. For cuprates, 
the model that serves as a sufficiently general microscopic starting point is the three band 
extended Hubbard mo deiiii. Each (planar) unit cell contains one copper 3d orbital and 
two oxygen 2p orbitals. The non-bonding combination of the 2p orbitals is not expected 
to play an important role. When this non-bonding combination is ignored, the three band 
extended Hubbard model reduces to a two band Anderson lattice like model. We will call 
this class of models as the two band extended Hubbard modelB. 

Inspired by these considerations, the theoretical question we ask is: can metallic non- 
Fermi liquids occur in the two band extended Hubbard model? We will address this question 
by incorporating general local interactions allowed by symmetry. Our goal is to treat inter- 
actions in a non-perturbative fashion, and seek to classify all the possible universality classes 
of this model. 

Recent work on this subject can be naturally separated into two categories. Work in 
the first category concern exclusively the single impurity physics. We have generalized the 
standard Anderson model by including all the on-site interactions allowed by symmetrylll^. 
The authors of Refs.ii"ii have generalized the standard Anderson model by introducing 
additional species of screening fermions. Work in the second category deal with the large D 
limit of the lattice extended Hubbard mo del&l^'S'i. Here we address the physics of the 
lattice model based on our understandings of the corresponding impurity problem. 

Connecting the local physics of an impurity model and a lattice model has a long tradi- 
tion. For instance, the slave boson large-A^ approach was first constructed to describe the 
Fermi liquid state of the single impurity Anderson modeli'i. The understanding of the 
impurity problem set the stage for the slave boson large- description of the coherent Fermi 
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liquid state of the Anderson-lattice modelU'El. 

B. Generalized Anderson model 

The generalized Anderson model we introduced0^0 is 

H = E^Ud + Uud^Udi + ^ Ekc\„Cka + 5Z ^{dlc„ + h.c.) 

ka f 

+Vndnc + JSd- Sc (1) 

Here, Uda = d^d^, Ud = J^a^da, = Ea4ca, 5"^ = (1/2) J2a,a' dlr^a' da', with r^., Ty, and 
Tz being the Pauh matrices, and Sc = (1/2) X^o-.o-' cJ-To-o-'Co-'. The first four terms describe 
the standard Anderson model. For the single impurity spin— | rf— orbital, the energy level 
is E^, and the on-site Coulomb repulsion U. For the most part, we will consider only 
the U = oo limit. For the spin— | conduction c— electrons, the energy dispersion is E^. 
c\ = ^j^^.^ J2k Cfco- is the Wannier orbital of the c— electrons at the impurity site. It hybridizes 
with the (i— electron through the hybridization matrix t. 

The last two terms are additional interaction terms allowed by symmetry. V describes 
a local density-density interaction between the impurity d— and local c— electrons. In the 
heavy fermion literature, this term is called the Falicov-Kimball interaction^. J describes 
the spin exchange interaction between the d— and c— electrons. It describes the sum of the 
direct exchange interaction and the indirect exchange interactions mediated by those high 
energy configurations not included in the model Hamiltonian. 

The standard Anderson model with a featureless conduction electron density of states has 
already been solved. In the particle-hole (p-h) symmetric case, namely for U + 2Ed = 0, the 
impurity c?— levels are illustrated in Fig. 3(a). For sufficiently large U, the empty impurity 
configuration (|0 >) and the doubly occupied impurity configuration (|2 >= | ti>) can be 
eliminated through a Schrieffer- Wolff transformatiorrll. The result is the Kondo problem 
with antiferromagnetic exchange interaction. The latter problem is so ved by a variety of 
methods, including scaling!2il numerical renormalization group (NRG)i, Bethe Ansatzi, 
slave boson large N methodE3, and conformal field theoryS. The conclusion is that, the low 
lying excitations can be well described by the strong coupling Fermi liquid fixed point. The 
local moment is quenched by the conduction electron spin polarization and hence disappears 
from the low lying excitation spectrum. 

In the p-h asymmetric case, U + 2E^ » \E^\, three impurity configurations have to 
be retained at low energies. This is the mixed valence problem. It differs from the Kondo 
problem in that low lying local charge (valence) fluctuations coexist with spin fluctuations. 
Historically, a variational study by Varma and YafetS, and RG studies of Haldane0 and 
Krishnamurthy et a/.l^, the Bethe Ansatz solutionis and the slave boson large-N resultsl^ 
have all found that the low energy behavior of the mixed valence problem is described by a 
strong coupling, Fermi liquid flxed point. This flxed point is qualitatively similar to that of 
the Kondo problem, though quantities such as the Wilson ratio are modifled. 

In Refs. we studied the generalized Anderson impurity model by extending Hal- 

dane's RG scheme such that the local charge fluctuations and local spin fluctuations are 
treated on an equal footing. In the mixed valence regime, there exist three, and only three. 
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kinds of fixed points. In addition to tfie aforementioned strong coupling Fermi liquid phase, 
there are two non-Fermi liquid phases which we termed the weak coupling phase and the 
intermediate phase. The strong coupling and weak coupling phases are the direct analog of 
the strong coupling phase of the antiferromagnetic Kondo problem and the weak coupling 
phase of the ferromagnetic Kondo problem, respectively. As for the local moment case, 
the weak coupling phase of the mixed valence problem requires that the exchange coupling 
be ferromagnetic. Therefore, this state is likely be of only academic value for the most 
part. The possible exception is the double-exchange model for the Perovskite Manganese 
OxidesS. The intermediate phase is unique to the mixed valence regime. Its existence 
came as a surprise. In this new phase, spin and charge excitations are separated; the spin 
susceptibility remains to have the Fermi liquid form as in the strong coupling phase, while 
the charge susceptibility and the single particle Green's function have an algebraic behavior 
with interaction-dependent exponents. Our RG results are substantiated by the strong cou- 
pling atomic analysis^ and by the exact solutions at certain exactly soluble points (Toulouse 
points)^. 

The single impurity model that Perakis et al.ii studied using NRG is defined as follows, 
H = E^Ud + Uud^ridi + Ekc\^Cka + Yl ^(4c<7 + h.c.) 

N N 

+Vndnc + J2 ^i^dnd + YY1 ^kclka^iM 

1=1 1=1 ka 

where cji^^, for I = 1,...,N, describe fermionic bands that interact, but do not hybridize, 
with the impurity ci— electron. The screening interactions Vi are introduced so that the 
Friedel sum rule is satisfiedi^. Ref. ^ reported NRG results in the mixed valence regime 
(with U = oo). The numerical results for the case of sufficiently large values of the screening 
interactions were interpreted as displaying divergent charge and spin susceptibilities near 
the mixed valence point. Such a phase is not expected from the Coulomb gas RG analysis. 
Within the Coulomb gas RGpicture, the effect of screening fermions is to modify the initial 
conditions of the RG The additional screening channels, while increasing the 

orthogonality effects, do not participate in the formation of fixed points other than those 
within our RG classification. For sufficiently strong V/, the Coulomb gas analysis predicts 
that the mixed valence state is the intermediate phase, with divergent charge susceptibility 
but regular spin susceptibility. While further NRG studies are clearly called for, here we 
note that the existing numerical data of Ref. might actually not be inconsistent with our 
Coulomb gas RG prediction. The reason is simple. Unlike for the charge susceptibility, the 
spin susceptibility continues to increase when is decreased through the transition regime. 
A true divergence is therefore hard to detect numerically. 
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C. Extended Hubbard model 



The two-band extended Hubbard model is the lattice analog of the generalized Anderson 
model Eq. (|ID, and is defined by the following Hamiltonian, 

H = Y1 ^d^di + ndi\ndil + Y ^iAaCja + Y K4aCia + h.C.) 

i i ij^cF ia 



+ ^{Vn-dinci + JSdi- Sci) (3) 

i 

The notations are essentially the same as in Eq. (|l|). The only difference is that we have 
used to label the electron level to emphasize the difference of this quantity in the lattice 
model with that of the impurity model (see the discussions around Eq. (|32D below). 

The first four terms describe the standard Anderson lattice model. The spin-| d— 
electrons have an infinite on-site Coulomb repulsion U, and an energy level e^- tij describes 
the kinetic energy term of the c— electrons. At every site, the d— and c— electrons hybridize 
with each other through the hybridization matrix t. The last two terms represent the on-site 
density-density and exchange interactions between the two bands at every site. 

Unlike for the single impurity problem, most of the studies on the Anderson lattice model 
in literature focus on the p-h asymmetric case. This is because the conventional heavy 
fermion metallic states are formed only in the p-h asymmetric case (the p-h symmetric 
case has received renewed interests due to the new developments on the so-called Kondo 
insulators^). We do not know as much about the standard Anderson lattice model as about 
the standard single impurity Anderson model. Only few of the theoretical methods that have 
been used in the single impurity problem are generalizable to the lattice case. This includes 
the Gutzwiller-like variational wavefunction£2lllll and the slave boson large-N methodSS. 
All these approaches have lead to the conclusion that the low energy regime is described by 
a Fermi liquid with heavy mass quasiparticles. The slave boson large N method has also 
been applied to the copper oxide modelB. In the metallic states without long range order, 
the solution is again a Fermi liquid. 

Our new understandings of the impurity physics, combined with the large D approach, 
have led to the conclusions that metallic non-Fermi liquid solutions are possible in the 
extended Hubbard model. This conclusion is firmly established in the large D limitEfflJii. 
In fact, it turns out that the mixed valence condition is much easier to realize in the lattice 
models than in the single impurity model. 



II. NON-FERMI LIQUIDS IN THE GENERALIZED ANDERSON MODEL: 

RENORMALIZATION GROUP 

We focus first on the single impurity problem. Given that the on-site repulsion is taken 
to be infinity, the generalized Anderson model can be thought of as a three level system 
with a particular form of symmetry breaking. The schematic picture is given in Fig. 4. The 
three levels correspond to the three impurity configurations: \a >= |0 >,| T>)| i>- The 
hybridization t, density-density interaction V, and spin-exchange J couple these three levels 
to the free conduction electron bath. Among the three levels, | t> and | |> are degenerate 
in the absence of external magnetic field. No symmetry, however, dictates the degeneracy of 
|0 > with I |>, I |>. It is therefore necessary to keep track of the energy level difference, E^, 
between | t>, | i> and |0 >. This section summarizes the RG analysis on the three- level 
system near its criticahty. Analysis of certain exactly-soluble points - the Toulouse points 
- is the topic of the next section. 
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A. Coulomb gas representation and renormalization group 



The RG analysis is based on a Coulomb gas representation0@ of the three-level system. 
This is carried out through an expansion in terms of the hopping amplitudes between the 
three configurations of the impurity problem. This is an extension of the classic work of 
Haldaneil such that the local charge fluctuations and spin fluctuations are treated on an 
equal footing. 

In practice, it is convenient to break the exchange coupling JSd ■ Sc into JzS^s^ + 
{J±/2){S^s~ + Sjs'^) where and J± represent the longitudinal and spin-flip components, 
respectively. V and are diagonal in the impurity configuration basis. Their effects are 
to cause different scattering potentials for the conduction electrons when they see different 
impurity configurations. When the impurity configuration is frozen in \a >, the scattering 
potential that the conduction electron of spin a experiences is V^: 

v:= V + j-jA 
v:= V - j,/4 

V,^= (4) 

Quite differently, the effects of t and J_l terms are to cause quantum transitions between 
different impurity configurations. Specifically, the hybridization t term causes transitions 
between the empty and singly occupied impurity configurations, and the spin-flip J±_ term 
the spin up and spin down impurity configurations. 

To construct the Coulomb gas representation, we expand the partition function in terms 
of t and Jx, and integrate out the conduction electron degrees of freedom. The resulting 
form for the partition function is a summation over histories, 

— = exp(-5[history]) (5) 

history 

where is the partition function of the free conduction electron sea. A history corresponds 
to a sequence of quantum mechanical hopping from one impurity state to another along the 
imaginary time axis. A transition from impurity state \a > to |/3 > is called a kink (a, (5). A 
history can be specified using the notation [ai, a„; ti, x„] which specifies an (a^, Oj+i) 
kink at time Tj. Fig. 5 illustrates a particular history. The statistical weight for a given 
history turned out to be 



Co 



+ Y^[K{ai, aj) + K{ai+i, a^+i) - K{ai, - K{ai+i, aj)]ln^— — - (6) 

^0 



i<j 



where ~ Po is the ultraviolet inverse energy cutoff. 

This action has the form of a Coulomb gas with two distinctive species of "Coulomb 
charges" . The two "Coulomb charges" correspond to the charge kink and spin kink, respec- 
tively. The fugacities of the two "Coulomb charges" are. 
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yt = yo,a = t^o 

Vj = Vhl = (7) 

As illustrated in Fig. 6, the charge fugacity yt corresponds to the hopping amplitude between 
two local states with different charge quantum numbers. Lilcewise, the spin fugacity yj 
describes the hopping amplitude between the | t> and | |> configurations. The fields ha 
describe the energy splittings among the three configurations. In the absence of an external 
magnetic field, ho = —^E^C,o and h^ = l-E'^^o- 

The logarithmic interactions between the hopping events reflect the reaction of the elec- 
tron bath towards the changes of the impurity configurations. The interaction strength is 
characterized by the stiffness constants, et = —K{0, cr) and ej = —K{'\, |) which in turn are 
determined by the bare interaction strength of the original Hamiltonian. Specifically, 

1 A<^ _ AO" _ 

Z 71 71 

e, = (1 - '-^-^f (8) 

7X 

where 5^ = tan~^(7rpoV^) is the scattering phase shift that the conduction electron bath of 
spin cr -whose density of states is po ^ experiences when the impurity configuration is frozen 
in \a >. 

The RG equations describe how the fugacities, stiffness constants, and the symmetry 
breaking field flow as we increase the cutoff ^. We follow the formalism of CardyEi. A 
detailed derivation can be found in Refs. Here we quote the results, 

dyt/dlni = (1 - et)yt + ytVj 
dyj/dln^ = (1 - ej)yj + yl 
det/dlni = -Qetyj + ejijj^ - y"^) 
dej/dlni = -2e,{yl + 2y]) 

dEai/dlni = {yj - y|) + Ed^l - 3y^) (9) 

The RG flow of the fugacities determine how the amplitude for making transitions between 
different impurity configurations are modified as we go towards longer time scales. When the 
amplitude grows the system is a Fermi liquid in analogy to the formation of Fermi liquid in 
the usual Kondo problem. When this amplitude renormalizes to zero, quantum coherence is 
destroyed and Fermi liquid theory breaks down. This way, we cast the breakdown of Fermi 
liquid theory in the framework of the macroscopic quantum coherence (MQC) problemB. 
The transitions between Fermi liquid and non-Fermi liquid phases are extensions of the well 
known localization transitions in the MQC problem with one essential difference. Here we 
deal with a special three-level system instead of the canonical two-level system studied in 
the MQC literature. This leads to a richer phase diagram that we describe below. 



B. Universality classes 

Solving the RG flow establishes the existence of three, and only three, mixed valence 
fixed points. The phase diagram for the mixed valence regime is specified in terms of the 
stiffness constants and ej and is given in Fig. 7. 
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The physical meaning of the strong couphng, weak couphng, and intermediate phases 
is most transparent when the three-level system is thought of as the hybrid of a two-level 
spin Kondo problem and a two-level charge Kondo problem. In the spin-Kondo problem, 
the "Coulomb charge" of the corresponding Coulomb gas corresponds to the spin-kink. A 
history of n spin kinks is shown in Fig. 8(a). The parameters of the Coulomb gas are the 
spin fugacity yj = Jx^nand the spin stiffness constant e'j = [1 — (2/vr) tan~^(7r^o<^2/4)]^. The 
RG flow is well known&i and is given in Fig. 9(a). For antiferromagnetic Jz, i.e. ej < 1, the 
flow is towards the strong coupling, Fermi liquid fixed point. While for ferromagnetic J^, 
i.e. ej > 1, the flow is towards a line of weak coupling fixed points. In a weak coupling fixed 
point, there is an asymptotically decoupled spin. 

The charge Kondo effect describes the physics of the so-called resonant-level model in 
the presence of a p-h symmetry. The spinless version of the Hamiltonian (|1|), with = 0, 
reduces to the resonant level model. It was realized some time agoElB that the resonant- 
level model can be asymptotically mapped onto the anisotropic Kondo problem, with the 
hybridization t and the density-density interaction V playing the role of J± and J^, respec- 
tively. The "Coulomb charge" of the corresponding Coulomb gas describes a charge kink. 
A history of charge kinks is illustrated in Fig. 8(b). The charge fugacity is yt = tC,o and 
the charge stiffness constant is = (1/2)[1 - (2/7r) tan-i«o1^/2)]^. The RG flow IS given 
in Fig. 9(b). The flow is towards a strong coupling Fermi liquid fixed point for < 1, 
but towards a line of weak coupling fixed points for > 1. In a weak coupling fixed point, 
the charge degree of freedom is asymptotically decoupled. Note that e[ > 1 corresponds to 
-V > -1/^"* = {2/7ipo)tan[{V2 - l)7r/2], i.e., a range of finite attractive density- density 
interaction. It is the charge analog of the ferromagnetic interaction. 

With this background on the charge sector alone and the spin sector alone, the meaning 
of the strong coupling phase of the mixed valence problem is transparent. Both the spin and 
charge Kondo problems are in the strong coupling regime; rapid fluctuations between all 
three local configurations take place and the conduction electrons quench both the charge 
and spin degrees of freedom of the impurity. It is expected that this phase is a Fermi liquid. 
A large-N analysis of this regime indeed gives rise to thiJli. So do the exact solutions of 
the two strong coupling Toulouse points (see next section). Likewise, in the weak coupling 
phase, neither the local charge nor the local spin degrees of freedom is quenched. Both 
the spin and charge Kondo problems are in the weak coupling regime, and all three atomic 
configurations decouple asymptotically at low energies. The weak coupling phase requires 
that the spin-exchange interaction be ferromagnetic. It is therefore very likely to be of only 
academic value with the possible exception of the double-exchange modefi. 

The unexpected phase is the intermediate phase. Here, the local spin degrees of freedom 
is quenched, but the local charge degrees of freedom is not. The charge Kondo problem 
is in the weak coupling regime despite of the fact that the spin Kondo problem is in the 
strong coupling regime. There are two local configurations carrying different charges which 
are decoupled asymptotically. The RG analysis establishes this as an allowed phase. The 
alternative situation, with the spin Kondo being in the weak coupling situation and at the 
same time the charge Kondo being in the strong coupling case, is not allowed: a relevant 
charge flipping (hybridization) drives the spin flipping relevant. The domain of attraction 
of the intermediate phase spans the parameter range within > 1 and ej < 1. The RG 
analysis cannot specify the precise boundary between the strong coupling and intermediate 
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phases; the dashed hne is only schematic. In terms of the parameters of the Hamiltonian 
(HI), with an on-site V term, this domain corresponds to a region with antiferromagnetic 
exchange interaction J and attractive density-density interaction V. Taking as an effec- 
tive parameter this condition can be satisfied in a variety of reahstic modelsE3. Finite range 
interactions also help realize these phases, as discussed in the context of impurity modelg^^H 
and in lattice model^i(see Section V). The transition between the different regimes is analo- 
gous to the localization phase transition studied in the context of the macroscopic quantum 
coherence problem0 and more recently in the context of transport through constrictions in 
interacting quantum wireJ^. 

There is one important difference in symmetry between the asjTiimetric Anderson model 
and the p-h symmetric resonant level model. The p-h symmetric resonant level model has 
a U(l) symmetry which ensures that the impurity level stays at the Fermi energy of the 
conduction electron sea. Equivalently, the singly occupied configuration and the empty 
configuration is guaranteed to be degenerate by symmetry. This degeneracy is responsible 
for the charge Kondo effect. In the case of the asymmetric Anderson model, no symmetry 
protects the degeneracy of the singly occupied and the empty impurity configurations. De- 
generacy can be achieved only through fine-tuning the bare impurity level. The condition 
for this degeneracy is none other than the condition for mixed valency. 

The phase diagram given in Fig. 7 applies only to the mixed valence regime. When the 
mixed valence condition is not satisfied, the impurity level is either too far below the Fermi 
energy or too far above the Fermi energy. They correspond to the local moment and empty 
orbital regimes, respectively. How far is too far away from the Fermi energy depends on 
whether the corresponding mixed valence state falls in the strong coupling, weak coupling, 
or intermediate domain. This crossover from local moment, mixed valence, to the empty 
orbital regimes are specified in the temperature versus impurity level space in Fig. 10. 

Fig. 10(a) specifies the finite temperature crossover for the strong coupling case. At zero 
temperature, the mixed valence crossover extends over a scale of ~ A*, the renormalized 
resonance width. The value of A* depends, of course, on where we are in the phase diagram 
Fig. 7. It is finite within the strong coupling domain. As we approach the phase boundary 
to the intermediate or weak coupling phases. A* vanishes in a Kosterlitz-Thouless fashion, 



A* ~ (po)~^exp[-l/Ve^"* - e] (10) 

where (e'^'^** — g) measures the distance from the phase boundary. 

For the intermediate phase. A* = 0, and the mixed valence point is a zero temperature 
critical point. At finite temperatures, there are three energy parameters: temperature (T), 
the running symmetry breaking field SEd{T) ~ {E^ - E^) - Ao(T^o)^^'*"^\ and the runnmg 
resonance width A(T) ~ Ao(T,^o)*"^'^*~^'' (where Aq ~ T^pot^ is the bare resonance width). 
The critical behavior associated with the mixed valence critical point occurs when \6Ed{T) \ < 
A(T) < T. This condition specifies the following crossover scale, 

1 



T'--m-E:i)/A,r^-^ (n) 

The correlation functions assume the form characteristic of the intermediate phase aXT > T' 
for a given E^, or equivalently, for a given temperature, when is tuned to the range 
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\E', - E2\ < Aoim)^'<-'^ (12) 

The intermediate mixed valence regime is a manifestation of the quantum critical phe- 
nomenon in the context of quantum impurity modelsEs. 



C. Further remarks 

Our most interesting finding in the mixed valence regime is the existence of a new 
phase, the intermediate phase. The Coulomb gas RG analysis provides the qualitative 
physical picture of the intermediate phase: spin excitations are quasiparticle-like, and charge 
excitations incoherent. However, the Coulomb gas RG analysis is not capable of determining 
the precise forms of the correlation functions for the intermediate phase (neither for the 
strong coupling phase, for this matter). This calls for alternative means through which 
correlation functions can be calculated explicitly. In the next section, we present explicit 
results of the correlation functions near several exactly soluble points. 

The Coulomb gas representation is based on the dilute instanton expansion. The RG 
analysis, while non-perturbative in the stiffness constants, is perturbative in terms of the 
fugacities. It is in principle possible that additional fixed points, not captured by the dilute 
instanton expansion, may occur. An example for the latter arises in the related, though 
qualitatively different, problem of tunneling through a point contact in a Luttinger liquidS'S. 
One way to probe the nature of the fixed points is to carry out a strong coupling atomic 
analysis, in the same spirit that Nozieres did for the usual Kondo problemE3. Namely, one 
analyzes whether the couplings are stable around the point where the couplings associated 
with the relevant fugacities take infinite values. We found that these points are indeed 
stable@, giving some support that the Coulomb gas RG classification of the universality 
classes are complete. This analysis of course does not completely rule out the existence of 
more fixed points. This issue became even more urgent due to the bosonization work of 
Refs. which reported fixed points unexpected from the Coulomb gas RG picture. It 



turned out that, as discussed in some detail in the next section, there are some technical 
subtleties with the application of the bosonization method to the mixed valence problem. 
When these subtleties are taken care of, the bosonization results become consistent with the 
Coulomb gas RG predictions. 



III. NON-FERMI LIQUIDS IN THE GENERALIZED ANDERSON MODEL: 

TOULOUSE POINTS 



There are three particular combinations of the interactiongl3 where the model is exactly 
soluble§. These points in the interaction parameter space are the mixed valence counterparts 
of the usual Toulouse point of the Kondo problemSS, and are naturally called the Toulouse 
points of the mixed valence problem. 

We identify the possible Toulouse points using the bosonization methoc@. Given that 
the interaction occurs at f = only, we need to keep only the 5— wave component of the 
conduction electrons. This S'— wave component is defined on the radial axis, r G [0,+oo), 
and can be further decomposed into an outgoing and an incoming components. In a standard 
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fashion, we extend to the full axis, x G (—00, +00), by retaining only one chiral component, 
which we denote by ipui^)- We can then introduce a boson representation for the ipui^) 
field. At the origin, 

^iix = 0) = Fj^e^*^ (13) 

Here, a is a cutoff scale which can be taken as a lattice spacing. $0- is the shorthand notation 
for $^(x = 0), 



= E + ^b,^e-^^/') (14) 




where bqa and are the Tomonaga bosons, and L is the length of the dimension and is 
taken to be infinite in the end. An important point is that, $0- depends only on the q ^ 
components of the Tomonaga bosons. In Eq. (|T^, the operator F^, and its adjoint F^-, are 
the so-called Klein factors. They should be thought of as acting on the q = sector of the 
Hilbert space for the Tomonaga bosons. More precisely, the Klein factors can be defined 
as the raising and lowering operators, respectively, in such a Hilbert spacdll'llB These 
operators are unitary, and anticommute among the different spin species. Furthermore, 
they commute with bg^ and 6^^. for g 7^ and, hence, also with <l>o-. 

The generalized Anderson model can be rewritten in the bosonized form, 

H = Hq + E'^'^ Xo-o- + H±t + H±j + Hv 

Hq = —[ rfx[(— )2 + {—f] 
An J dx dx 

V 27ra o- 

where we have used Ud = T,aXaa, dl = X^q, S'j" = X^^, and = (X|| - X^)/2. X^/s = 
\a >< f3\ are the Hubbard operators. The requirement that a,l3 take three, and only three, 
impurity configurations, |0 > and \a >= d}^\0 >, amounts to the following constraint, 

Xtt+Xu+Xoo = l (16) 



In Eq. ([T5|) , $c,s = i^^^^i)/ are the charge and spin bosons, respectively. = ^ J2a 
and = 75 ~ l/27rpo is the Fermi velocity. 

The Toulouse points are derived through applying a canonical transformation to the 
Hamiltonian Eq. (|T5|) and demanding that the transformed H±t and H±j have simple forms 
and, simultaneously, the transformed Hy vanishes. Three such Toulouse points exist. The 
details are given in Ref. In the following, we only quote the effective Hamiltonian, and 
the results for the single particle, spin-spin, and charge-charge correlation functions, at each 
of these Toulouse points. 
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A. Strong coupling Toulouse point I 



The first Toulouse point corresponds to ej = and ej = 0. According to the phase 
diagram (Fig. 7) of the Coulomb gas RG analysis, this point lies deep in the strong coupling, 
Fermi liquid region. 

To write the effective Hamiltonian, we need to introduce pseudoboson operators fej. and 
bl defined as follows, 

Xaa' = ft fa' 

Xao = f% 
Xqo = blbo 

bi = fM (17) 

Note that the pseudoboson operator incorporates a Klein operator associated with the 
conduction electron degrees of freedom. In terms of these operators, the constraint Eq. (|TBp 
can be rewritten as X^o- b\ba + &o&o = 1- The effective Hamiltonian can then be conveniently 
written as 

H,, = E°(^ bX - bib,) + Y.(b% + H.c.) - ^{b\b^ + H.c.) 

a- V27ra „ 47ra 

^H={^){Y.b%-hlb,){Ui^'^- 



ZTvpo ^ ZTT ax 

(18) 

where Kc and measure the deviation from the Toulouse point. This effective Hamiltonian 
is composed of three parts: H^i is the Hamiltonian for the isolated three levels, ac >, 
ac >, and bllvac >, where \vac > denotes the vacuum state. The t and J± are transverse 
fields, and provides a longitudinal field. Hq describes the free spin and charge bosonic 
fields. Finally, AH is the dissipative term coupling the three levels to the bosonic bathes. 
The effective Hamiltonian therefore is a three- level generalization of the two- level "spin"- 
boson problemil0. 

All the correlation functions of this three- level "spin" -boson problem can be calculated 
explicitly. The results for the single-particle Green's function Gdir) = — < Trda{T)dl{0) >, 
the density-density correlation function Xp(t) =< TT-nrf(r)nrf(0) >, the longitudinal and 
transverse spin-spin correlation functions Xct^(t) =< T^S'^(r)S'^(0) > and xt'i'^) 
Tt-S~{t)S^{0) > are given as follows, 

G.(r) 



r 
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where hg = J_|_/47ra and = t/ \/27ia. A long time 1/r behavior for the single particle 
Green's function, together with the behavior for the two particle correlation functions, 
imply that the system is a Fermi liquid. Therefore, this Toulouse point describes the strong 
coupling phase. 

Unlike for the Kondo problem, keeping track of the anticommutation relation between 
fermions of different spins in the boson representation plays an essential role. Had we not 
included the Klein operator in the boson representation of the fermion operator Eq. (|13|), 



the pseudoboson operator b^^ defined in Eq. (0) would not include the additional Klein 
operator (it would then be a pseudofermion, as a matter of fact). The sign of the J± term 
would be reversed. It can be seen by diagonalizing the three level atomic problem, H^i, 
that a level-crossing would arise as is varied. The critical value of where levels 
cross would correspond to a non-Fermi liquid critical point. A signature for the unphysical 
nature of the non-Fermi liquid critical point associated with the level crossing can be seen by 
comparing the transverse and longitudinal spin-spin correlation functions. It can be shown 
that the longitudinal spin-spin correlation function has a non-Fermi liquid form, but the 
transverse spin-spin correlation function retains the Fermi liquid form. As a matter of fact, 
the level crossing would have occurred had we started from an unphysical model with 
antiferromagnetic but J± ferromagnetic. 

Within the bosonization approach, the meaning of the atomic configurations in the canon- 
ically transformed bases is somewhat obscure. The physical content of these configurations 
becomes transparent once we compare them with the atomic configurations that appear in 
a perturbation expansion of the original Hamiltonian in terms of J±/Jz, J±/V, t/J^, t/V, 
W/Jzi and W/V. This atomic analysis is carried out in Ref. H3, from which it is physically 



clear that the ground state is always a singlet no matter what the value of E^ is. No level 
crossing is expected! 

B. Strong coupling Toulouse point II 

This corresponds to = 1/2 and ej = 0. The Coulomb gas analysis would again predict 
this point to be deep in the domain of attraction of the strong coupling Fermi liquid phase. 
The effective Hamiltonian is, 

H^ff = Ho + E", Y: flL + ^ [(E ~fl)^ + H.c] 

Here, rjl denotes a spinless conduction electron band; it comes from refermionizing the charge 
boson, = Xo-qF^FJ^ is a pseudofermion operator. Unlike for H^jj-, we have kept only the 
Kg term in AH, as the term is not important. All the correlation functions can once again 
be explicitly determined, 

r 

Xa M~ ( — ) 
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2TipQhs T 



(21) 



where Ks is again the deviation from the Toulouse point. Once again, the long time behavior 
of the single particle and two particle correlation functions establishes the strong coupling, 
Fermi liquid nature of this Toulouse point. 

Except for the change of sign in J±, the effective Hamiltonian (|20|) is identical to that 
of Refs. ^ and |5l|. Once again, when the Klein operators are properly incorporated in the 
bosonization representation of the fermion fields, no level crossing occurs. 



C. Toulouse point for the intermediate phase 

This last Toulouse point occurs at et = 1 and ej = 0. It is not inconsistent with the 
Coulomb gas results that these values of the Coulomb gas stiffnesses lie close to such a bound- 
ary (though it is not possible to determine the precise boundary between the intermediate 
phase and the strong coupling phase from the Coulomb gas analysis). 

In order to write down the effective Hamiltonian in a convenient fashion, we need to 
introduce a new basis set for the three levels, 1^4 >, \B >, and |0 >. They are defined as 
follows, 

\B> = ^j:i-flH)\vac> 

\0> = bl\vac> (22) 

where /j^ and bl are pseudofermion and pseudoboson operators defined in Eq. (pT]). In this 
new basis, = flfa = {Xaa + Xbb - Xqo), = (1/2) (^flfa = {Xab + Xba)/2, 
and X^iFfF^ = —Xaa + Xbb- The effective Hamiltonian has the following form, 

^eff = H ^kcl^^kcT + 2ty/na[XA0C^Ci + H.C.] 
fccr 



— {Xaa + Xbb - Xoo)(c|c| + cjq) 



2vrpo 



+^{Xab + Xba){c\c^ - cjq) (23) 

In this effective Hamiltonian, the charge sector is described by a genuine charge Kondo 
model. \A > and |0 > play the role of | t> and | |> of the anisotropic spin Kondo problem 
and should be thought of as objects carrying charge 2 and 0, respectively. The transformed 
hybridization term is the direct analog of the spin-flip term in the anisotropic spin Kondo 
problem. The residual interaction in the charge sector, is the analog of the longitudinal 
exchange interaction in the anisotropic spin Kondo problem, with the density playing the 
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role of the spin in the latter. The essential difference between the charge Kondo problem in 
this mixed valence context and the spin Kondo problem lies in the symmetry-breaking field. 
In the latter, the spin rotational invariance guarantees that no explicit magnetic field term 
will be generated in the absence of an external magnetic field. In our charge Kondo problem, 
the p-h symmetry is explicitly broken, and the symmetry-breaking field h^^^^^^ — E^ — ^ is 
in general non-zero. For the impurity problem, the condition that the renormalized h^^'^^sf: 
vanishes can be achieved only through fine-tuning the bare d— level to a critical value 
El 

When /i^harge — Q is enforced, a zero temperature quantum phase transition takes place 
increased through zero. The transition is characterized by a Kosterlitz-Thouless 
transition in the charge sector; the spin sector is not critical. The phenomenology of the 
intermediate phase is recovered on the negative Kc side, to which the remaining of this section 
is devoted. Here, the charge sector is described by the weak coupling fixed points of the 
charge-Kondo problem, while the spin excitations by the strong coupling, Fermi liquid-like 
fixed point of the Kondo problem. A spin-charge separation takes place. 

Within the charge sector, the impurity configuration in the ground state is entirely |0 > 
for h'^^^'^^'^ < 0, and \A > for h'^^^^'^^c ^ This is the result of infinite charge susceptibility 
in the corresponding ferromagnetic charge Kondo problem. Exactly at h'^^'^'^^e _ namely, 
when E^ is tuned to the critical value E^ = the impurity degrees of freedom in the 
ground state involve an equal, incoherent, mixture of |0 > and 1^4 >. Schematically, the 
ground state wavefunction can be written as0=|A>0^ + |O>0o where (pA and 0o are 
the wave functions of the conduction electrons such that is the solution to a ferromagnetic 
Kondo model with zero magnetic field. With h^^^'^s^ = 0, the intermediate mixed valence 
dynamics applies at all temperatures. When E^j^ is moved away from the critical value, a 
finite cross-over temperature T^o ~ \E^ — E^\ emerges. The intermediate mixed valence 
dynamics continue to apply at T > Tco- At low temperatures {T < Tco), however, the charge 
fluctuations become gapped out. 

The singlc-particlc, density-density, longitudinal and transverse spin-spin correlation 
functions are given as follows. 



X-{r)- {^n'-^f (24) 

where Xp labels the connected part. The exponent for the single particle Green's function 
is particularly noteworthy. The | part is the contribution of the spin degrees of freedom. It 
is independent of interactions, and is the same value as wc would get for a non-interacting 
problem! The remaining part, |(1 — V^^Y is due to the charge degrees of freedom and is 
interaction dependent. This is consistent with the physical picture that in the intermediate 
phase the low lying spin excitations are quasiparticle-like while charge excitations have the 
non-Fermi liquid form. 
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Other correlation functions in the charge sector also have an algebraic behavior with 
interaction-dependent exponents, and a pairing susceptibility, 

< T.(Ec.rf,)(r)(E44)(0) >~ (^)(^)' (25) 

is enhanced compared to the Fermi-liquid case. This makes it plausible that the ground 
state in the corresponding lattice model is superconducting. In that case, the intermediate 
mixed valence dynamics would describe the physics in the normal state, i.e., at temperatures 
between the transition temperature and some upper cutoff energy scale. 

To summarize, the explicit results for the correlation functions in this Toulouse point 
highlight all the features expected of an intermediate phase: spin-charge separation; a quasi- 
particle residue vanishing in a power-law fashion; Fermi liquid like spin correlation functions; 
and self-similar local charge correlation functions with interaction-dependent exponents. We 
note in passing that these characteristics bear strong, similarity to those of the Luttinger 
liquid in one dimensional interacting fermion systemJ^^^^^ 



IV. THE EXTENDED HUBBARD MODEL AS A LATTICE OF ANDERSON 

IMPURITIES: LARGE D LIMIT 

We now turn to the lattice model, Eq. (^), which we suggestively rewrite as 

i <ij> 

hi = e°^ndi + Undi^Udii + t(J2 4aCia + h.c.) 

cr 

+Vndina + JSdi ■ Sci 
hij = tij ^ cfg-Cjo- (26) 

cr 

This is pictorially illustrated in Fig. 11(a), in which each black dot represents an hi term, 
and across each bond on the lattice there is an hij term. 



A. Mapping to a self-consistent impurity problem in the large D limit 

The limit of infinite dimensionsilS is defined by scaling the hopping term, tij, in terms of 
the dimensionality {D) such that the limit is well-defined. For the nearest neighbor hopping 
term 

t<ij> = (^''') 

The D oo limit is taken with kept fixed. 

We recall that, the large D limit of a classical non-frustrated lattice spin system is 
taken by scaling the nearest neighbor coupling to be of order 1/D. For any given site, the 
homogeneous contributions from neighboring sites add up to an effective field, of order unity, 
that acts on the spin of the selected site. All other contributions are of finite orders in 1/D 
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and vanishes in the large D hmit. This is the content of the Weiss molecular field theory for 
classical magnets. In the quantum systems, the single particle hopping contributes to the 
kinetic energy of the fermion system which, at zero temperature, is the zero point energy 
associated with quantum fiuctuations. The 1/a/D scaling in Eq. (|27|) , as opposed to the 
1/D scaling, is necessary to capture these quantum fiuctuations. With this scaling, the 
average kinetic energy per unit cell is of order unity in the large D limit. 

In finite dimensions, when on-site interactions are present, a single partile hopping term 
will generate effective intersite interactions involving two or more particles. With the single 
particle hopping term being scaled as in Eq. (^), the generated interactions are of order 
1/D or higher. 

For a selected site, say site 0, the effect of the rest of the sites is to generate a retarded 
Wiess mean field that couples to the single particle degrees of freedom at site zero. The 
modifications to the on-site dynamics involving two or more particles are higher order in 
1/D and do not survive the large D limit. The result is that, all local correlation functions 
of the lattice model can be entirely determined by the following effective on-site action: 

Stli = So- dr dr' Y: cUr)g,\r - T')c,.{r') (28) 

5*0 is the action associated with h^. Since contains all the local interactions, the proce- 
dure treats the local interactions in a dynamical fashion, g^^ir — r'), or equivalently, its 
Fourier transform, gQ^{iuJn), where is the fermionic Matsubara frequency, is retarded. 
This is the result of integrating out the rest of the degrees of freedom other than site 0, 
at the one particle level. Pictorially, gQ^ describes the effect of all the Feynman trajecto- 
ries in which one electron leaves site zero, explores the lattice, and returns to the origin. 
Translational invariance demands that the local correlation functions of the lattice model 
are site- independent, and are the same as the correlation functions of the impurity model. 
This leads to the following self-consistency equation, 

g^^{iun) = -^^UotojlGijiiUn) - Gio{iuJn)Goj{iLJn) / Goo{iuj„)] (29) 



where Gim{T) = — < TrCia{T)cl^^{0) >H is the lattice Green's function. Eqs. (|28| , |29|) define 
the dynamical mean field formalism that is exact in the large D limitil. 

It is physically more transparent to rewrite S'^^'^ in the Hamiltonian form. We achieve this 
by introducing a non- interacting electron bath whose dispersion and coupling to the cqo- has 
to be determined self-consistently. Introducing r]l^ and r]ka as the creation and annihilation 
operators for this fermion bath, where A; is a dummy momentum variable, we can rewrite 
the effective impurity problem in terms of the following effective impurity Hamiltonian, 

H!li = ho + J2 tkiviaCoa + H.c.) + ^ ek7]l^r]ka (30) 

kcr ka 

The self-consistency equation (^) is equivalent to 

Ytl/{iujn - €k) = -^^tiotojlGijiiuJn) " Gio{iLJn)Goj {iuj„) / Goo{iuJn)] (31) 

k ij 

Pictorially, we have reduced the task of solving the full lattice problem of Fig. 11(a) into 
solving a fully interacting quantum impurity embedded in a self-consistent fermionic sea, as 
is illustrated in Fig. 11(b). 
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B. Solution to the effective impurity problem 



The crucial question is the nature of the self-consistent bath. What we found is that, 
as long as the solution is metallic, the density of states of the fermionic bath at the Fermi 
energy is finite. This is a self-consistent statement. The reasoning goes as follows. Assuming 
that the bath density of states is finite at Fermi energy, we can proceed to solve the impurity 
problem in an asymptotically exact fashion by applying bosonization technique. Among the 
quantities that can be calculated asymptotically exactly is the local Cq Green's function. 
That this Green's function has a regular spectral function is seen, in the Coulomb gas 
representation, by noting that the local Cq does not creat a kink. A regular cq Green's 
function implies a regular self-energy for the c— electrons, which, through the self-consistency 
equation (PT]), in turn implies that the density of states of the self-consistent fermionic 
bath is regular! The self-consistency is hence established. This argument applies to any 
lattice. In the special case of Lorentzian density of states, this statement is more than 
asymptotically exact; it is exact. In the case of a Bethe lattice with infinite coordination, we 
have numerically solved the self-consistency equation^. We indeed found that the density 
of states of the bath fermions at the Fermi energy is finite as long as the solution is metallic. 
In Fig. 12, we plot the imaginary part of the c?— electron and c— electron Green's function 
as a function of the Matsubara frequency. The zero frequency limit of this Green's function 
is identical to the density of states at the Fermi energy. It is clear that, even though the 
(i— electron Green's function is divergent, as expected in the non-Fermi liquid form discussed 
in the previous sections, the c— electron density of states is regular. 

Armed with this understanding of the fermionic bath, the only essential difference of 
this self-consistent Anderson model with the single impurity generalized Anderson model is 
that, this time the effective ci— level is 

Er = e',-H (32) 

instead of e°. The effective impurity level here is measured with respect to the Fermi energy 
of the lattice model which is, of course, different for different amount of electrons. 

The fact that the density of states of the self-consistent fermionic bath at the Fermi 
energy is finite implies that, we can carry out an asymptotically exact analysis on the 
self-consistent generalized Anderson impurity model exactly the way we did for the single 
impurity generalized Anderson model. The classification of the possible phases of the single 
impurity Anderson model applies to the self-consistent Anderson impurity model, provided 
that we express the stiffness parameters of the phase diagram. Fig. 7, in terms of the 
self-consistent parameters. In particular, there is a metallic non-Fermi liquid state of the 
extended Hubbard model that corresponds to the intermediate phase of the impurity model. 
As in the impurity model, we have a spin excitation spectrum that is spin-| quasiparticle- 
like, and a charge excitation spectrum that is incoherent. What we have is a local route 
towards spin-charge separation in the extended Hubbard model. 

The fact that the impurity level of the self-consistent impurity problem is measured with 
respect to the chemical potential implies that, for a model with fixed e[J, the temperature- 
chemical potential crossover is given in Fig. 13(a). 
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C. Pinning of chemical potential 



A remarkable phenomenon takes place. This is the pinning of chemical potentialE3. The 
critical chemical potential, at which the mixed-valence state persists to zero temperature, 
corresponds to a range of electron densities. It can be seen as follows. The local correlation 
functions of the extended Hubbard model in infinite dimensions are given by the impurity 
problem. In particular, the occupation numbers of the lattice, for a given chemical potential, 
can be obtained from the local Green's functions of the corresponding impurity model. It 
follows from our analysis of the impurity model that, at zero temperature, ri^ (and also 
n = nd + ric) are discontinuous functions of the chemical potential: ~ nj" ^ 1 — O(t^) for 
/i > fic, while rid ^ ^ < f^c- At finite temperatures, na is increased from 

to as /i is increased from fic — A/i to fic + A/i where 

A/i ~ AoiTpoY^'*^-'^ (33) 

As long as is over the range {n^,n^), the condition Eq. (p!2D is satisfied. The correlation 
functions will be controlled by the intermediate phase at criticality. 

It is remarkable that the fact that our impurity model is associated to a lattice problem 
forces the effective impurity model to be at criticality, with a larger symmetry than we would 
have naively expected. Physically, for an incoherent state to be metallic, it is necessary to 
allow charge transfer between the local degrees of freedom and the bath. This can only 
happen if the local charge degrees of freedom is in equilibrium with the conduction electron 
bath. This requires the heavy level to be at the chemical potential. 

D. Further remarks 

The RG analysis that leads to the classification of the possible phases of the effective 
impurity problem discussed in Section II is based on a small hybridization expansion. In 
the context of the low energy effective Hamiltonians for the high system, one of the 
important questions is whether the extended Hubbard model can be further reduced to an 
effective one band Hubbard modelSS. One argument used in this context is that, when 
the hybridization is sufficiently large, it is more appropriate to first diagonalize the problem 
within a unit cell, leading to the Zhang-Rice singletE3. The effective Hamiltonian for the 
low energy local orbitals is the so called t — J Hamiltonian. In light of this construction, 
a natural question to ask is whether the physics of the extended Hubbard model at large 
hybridization is different from that at small hybridization. The large D limit provides 
a unique opportunity to address this question. This problem has only been numerically 
studied in the spinless version of the extended Hubbard modeS. The numerical solution 
indicates that the qualitative phase diagram is similar for the large and small hybridization 
limits. However, the precise values of the exponents of the correlation functions could be 
modified as the hybridization is increased. Further work along this direction needs to be 
carried out. 

We have established the existence of metallic non-Fermi liquid states in the extended 
Hubbard model Eq. (^ in the large D limit. What happens in finite dimensions? This is a 
question which is only beginning to be addressed. Some progresses are reported in the next 
section. 
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We close this section on a methodological note. Various numerical methods, such as 
quantum Monte Carlo&^ and exact diagonalizationSS, can be used to solve the self- 
consistent dynamical mean field equations associated with the D = oo limit. Whatever the 
means, the solution to these equations should always describe the solution to an impurity 
coupled to a self-consistent fermionic bath. And it is always instructive to ask a) what is 
the nature of the density of states of the self-consistent fermionic bath near its Fermi energy 
(is it regular, gapped, vanishing with a power law, or singular with a power law, to name a 
few); and b) what are the low lying levels associated with the impurity. Armed with these 
information, it is usually possible to use RG or other analytical means to determine the 
qualitative behavior of the solution. 



V. COMPETITION BETWEEN THE LOCAL AND SHORT RANGE 
FLUCTUATIONS: ALTERNATIVE LARGE-D LIMIT 



One major advantage of the large D approach is being able to treat local correlations 
in a fully dynamical fashion. This feature is responsible for our uncovering the metallic 
non- Fermi liquids in the extended Hubbard model that other methods have failed. One ma- 
jor disadvantage of the large— D approach is that, spatial fluctuations beyond one particle 
level are all frozen: intersite interactions reduce to Hartree contributions only. For physical 
systems in finite dimensions, intersite RKKY or Superexchange type interactions are ex- 
pected to compete with local correlations. For instance, an unstable non-Fermi liquid fixed 
point arises due to the competition between the inter-impurity RKKY interaction and the 
local Kondo couplings in the two-impurity Kondo problem^. In the Kondo lattice models, 
such a competition led to the competition between long range magnetic ordering and Kondo 
singlet formation. In the absence of long range ordering, the dynamical role of the intersite 
interactions on the local Kondo-like physics has largely been left unexplored. From the 
large— D point of view, one way to recover the spatial fiuctuations is the perturbative 1/D 
expansion. Truncating the perturbation series to order requires solving at once clus- 

ters containing one, two, ...,n + 1 sites embedded in their respective self-consistent media@. 
The practicality of this procedure is unclear at this stage. An alternative route is a loop 
expansion with the requirement that the D = oo results be recovered at the saddle-point 
level, as has been constructed in models with certains forms of quenched disordert3. For 
clean systems, it turns out to be difficult to formulate such a loop expansion. 

We have recently introduced an alternative large D limit to study the interplay between 
local correlations and short range spatial fiuctuations in the two band extended Hubbard 
modelii'S. In this procedure, an explicit intersite density- density interaction term is intro- 
duced, and is scaled in terms of the dimensionality such that its fluctuation part survives the 
large D limit. This procedure leads to an impurity embedded in a self-consistent fermionic 
bath and a self- consistent bosonic bath. Detailed analysisll has so far been carried out only 
for the spinless version of the extended Hubbard model, given by the following Hamiltonian, 

H = Y^lE^Udi + t{d\ci + h.c.) + Vndinc,i\ 

i 

+ Yl l^ij'^Uj + % ■iT'di-- ndj :] (34) 

<ij> 
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The standard large D hmit is taken with tij of the form Eq. (pT]) and Vij of order 1/D. 
With that scahng, only the static component of Vij gives a non-vanishing contribution in the 
large D limit. Hence, intersite interactions give only a Hartree contribution. The alternative 
large D limit is taken with Eq. (pTl) and 

v<ij> = Vo/VD (35) 

In order to have a well-defined large D limit, it is necessary that the zero frequency mode, 
i.e. the Hartree term, be treated separately. In the absence of symmetry breaking, the effect 
of the Hartree term is a change to the chemical potential. This is handled through normal 
ordering, : n := n— < n >. 

The procedure outlined in the previous section leads to the following effective impurity 
action. 

In addition to the self-consistency equation (^), an additional self-consistency equation is 
required, this one for the density propagator, 

Xo^iii^n) = Y^VioVojixijiii^n) " Xioiii^n)Xojiii^n) / Xooiii^n)] (37) 

where Xim{i^n) is the Fourier transform of the lattice density- density correlation function, 
Xim{T) =< Tr : ni : (r) : : (0) >h- Wn is the bosonic Matsubara frequency. 

The effective action can once again be written in terms of a single impurity Hamiltonian. 
In our spinless case, this is a self-consistent resonant-level model with an additional screening 
bosonic bath, 

HtUp = (e° - 1^)^^ + tidlco + h.c.) + VudoUco 
+ J2 tkivlco + h.c.) + J2 ^kvlvk 

k k 

Here, rfl. creats, like in the previous section, a fermionic bath with a dummy momentum 
variable k. The dispersion, e^, and the hybridization coupling parameter, tk, are to be 
determined self-consistently. Likewise, creats a bosonic bath with a dummy momentum 
variable q. The corresponding self-consistent dispersion and coupling parameters are Wq 
and Fq. 

Detailed analysis shows that, in this case, the fermionic bath density of states remains 
regular. The spectral function of the bosonic bath, however, can be highly non-ohmic. This is 
fortunate, for an impurity model with a fermionic bath having an arbitrary form of density of 
states near the Fermi energy is quite difficult to handld^. On the other hand, an impurity 
model coupled to a bosonic bath with non-ohmic spectral function can still be analyzed 
asymptotically exactly, within a modified kink-gas picture. Details are given in Ref. 
The most interesting regime is again the mixed valence regime, for which the renormalized 
effective impurity level is zero. The self-consistent solution is shown in Fig. 14. The phase 



23 



diagram is specified in terms of tliree parameters, gt = pot, gv = [I — {2/7i)tan~^{7ipQV/2)], 
and g^ = PqVq. Tliey are essentially the dimensionless hybridization, on-site density-density 
interactions, and intersite density- density interactions. 

At Op = 0, the problem reduces to the usual large D case discussed in the previous 
sectionEj. For our spinless problem, a Kosterlitz-Thouless transition takes place describing 
the charge-Kondo effectJll. When gy < fi'y**, i.e., —V < Vq"^, the solution is a Fermi liquid. 
For gv > S'y**, i.e., —V > Vq"^, the solution is a line of non-Fermi liquids with the connected 
local density susceptibility. 



The exponent a is interaction-dependent, increasing from to 2 as one moves away from 
the critical pointc29'E21. 

The intersite interaction vq modifies the phase diagram in several ways. Consider first 
9v > S'y**- The line of fixed points of the f o = problem becomes unstable. Remarkably, 
the correlation functions in the new fixed points can be determined. In fact, they have the 
same form as given in Eq. (^). 

For gy < gy^^, we are able to establish the existence of a phase transition as Vq is 
increased. For sufficiently strong vq, the solution must be a non- Fermi liquid. Physically, 
the intersite density-density interactions provide charge-screening, which contribute to the 
orthogonality effect. In the mixed valence regime, this orthogonality helps realize the weak 
coupling fixed point with incoherent charge excitations. For sufficiently small Vq, on the 
other hand, the Fermi liquid solution is stable. 

As a result, non- Fermi liquids with self-similar correlation functions occur even for repul- 
sive values of the on-site density- density interaction. It is therefore not necessary to require 
attractive on-site interactions to realize the incoherent charge state. 

A finite intersite interaction vq also changes the nature of the phase transition. We have 
shown thatH the phase transition is not of the Kosterlitz-Thouless type. We have been 
unable to establish the precise nature of the phase transition. 

The self-consistent equations and the kink gas analysis can also be carried out in the 
spinful Hubbard model. The form of scaling equations implies that the results derived 
here for the spinless model carries over to the charge sector of the spin-charge separated 
intermediate phasell^i. 

The formalism outlined in this section can be generalized to various different contexts. 
The effects of intersite spin exchange interactions are the obvious next problem to study. 
Three or more particle intersite interactions can also be treated along this line. 

VI. CONCLUSIONS, NEW INSIGHTS, AND OPEN QUESTIONS 

The focus of this review article has been on the competition between local charge and 
local spin fluctuations, both in the single impurity Anderson model and in the lattice ex- 
tended Hubbard model. We have found that such a competition leads to metallic non-Fermi 
liquids in certain interaction parameter range. In particular, we have identified a novel 
non-Fermi liquid mixed valence state, called intermediate phase. This phase displays the 
phenomenon of spin-charge separation. The low energy spin excitations are much like that 




(39) 
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of the strong coupling Fermi liquid phase as one would derive from, for instance, the slave 
boson condensed phase within the slave boson large-N approach. The charge excitations are 
distinctively of a non-Fermi liquid form. 

Are these non-Fermi liquid states relevant to real materials? For single impurity prob- 
lems, these non- Fermi liquid phases can be realized only when the impurity level is tuned to 
be close to the Fermi energy of the conduction electrons. We can envisage two contexts in 
which this kind of fine tuning of the impurity level is physically feasible. The first is in the 
context of dilute impurities in metals. The Fermi level of the conduction electron sea can be 
varied by substituting some of the elements in the compound with ones of different valency. 
Called Fermi level tuning, this mechanism has already been invoked to explain the trend of 



the Kondo energy in certain Uranium based heavy fermionalH^'tHJ. Valence fluctuations in 
Uranium based compounds are in general much stronger than in Cerium based compounds. 
It is conceivable that some of them have interaction parameters that fall in the domain of 
the intermediate phase. Systematic studies of the Fermi level tuning phenomenon, therefore, 
can play a significant role in the current search and study of non- Fermi liquids in /—electron 
based materials. The second context where impurity level can be tuned is in mesoscopic 
systems. This time, the tuning is achieved through biasing the confined area with respect 
to the leadsEsi. 

For lattice problems, the level tuning requirement is much less stringent. This is the 
result of the phenomenon of the pinning of chemical potential, discussed extensively in Sec. 
IV. There is a range of electron density over which the effective impurity level lies close 
to the Fermi level of the self-consistent fermion bath. Our findings of the existence of 
both the strong coupling and intermediate phases in the mixed valence regime provide new 
insights into the similarities and differences between the heavy fermions and high cuprates, 
mentioned at the beginning of this manuscript. The spin dynamics of the intermediate 
phase and the strong coupling phase are similar, both displaying a crossover from the high 
temperature local moment regime into the low temperature coherent regime. The charge 
dynamics, on the other hand, are qualitatively different in these two phases. In the strong 
coupling Fermi liquid phase, the charge dynamics track with the spin dynamics. In the 
intermediate phase, the charge dynamics have non-Fermi liquid behavior characterized by 
the correlation functions discussed in Sections II and III. 

The phenomenology of the conventional heavy fermions, such as CeCug and UPts, are 
well described by the strong coupling Fermi liquid phase. Can the normal state of the high 
Tc cuprates be described by the intermediate phase? The contrasting behaviors seen in the 
temperature dependences of the NMR relaxation rate and the electrical resistivity in the 
cuprates (Fig. 2), when viewed in the context of those of the heavy fermions (Fig. 1), are 
consistent with the quahtative differences between the spin and charge dynamics expected 
in the intermediate phase. However, these two quantities measure very different correlation 
functions. The NMR relaxation rate measures mainly the momentum {q) averaged electron- 
spin response, while the electrical resistivity the q = electrical current-current correlation 
function. Therefore, the precise implications of the contrasting temperature dependences 
of these two quantities in the cuprates are hard to specify. One clear-cut probe of the 
relationship between the spin and charge excitations would be to compare the temperature 



dependences of the electron-spin resistivity and the electrical resistivityll^. This requires an 
experimental measurement of the electron-spin diffusion constant. 
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Our results also raise a number of theoretical questions. Unlike for the multi-channel 
Kondo problem, the lattice model we have studied, Eq. (|^), has a well defined limit of 
vanishing interactions, U ^ 0, V ^ and J ^ 0. In this limit of vanishing interactions, and 
for dimensions higher than one, the perturbative RG analysisEl would identify no instability 
towards a metallic non- Fermi liquid state. By focusing on the strong coupling limit, U = oo, 
and taking the limit of infinite dimensions, D = oo, we have identified non- Fermi liquid 
solutions. Obvious questions arise: a) what happens as the on-site Hubbard interaction U 
gradually decreases from infinity, all the way to f/ = 0? If a phase transition takes place, is 
it also of the Kosterlitz-Thouless form? b) What happens when the dimensionality decreases 
from infinity to physical dimensions? The approach outlined in Section V provides a starting 
point to address one aspect of the finite dimensionality effects, namely the competition 
between the on-site and inter-site correlations. The results summarized there imply that the 
non-Fermi liquid phases survive the short range spatial fluctuations. The critical behavior 
of the quantum phase transition from the Fermi liquid to non-Fermi liquid states, on the 
other hand, are strongly modified by the spatial correlations. As for a), it remains an open 
problem at the present time. 
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FIGURES 




FIG. 1. The NMR relaxation rate (1/Ti) and electrical resistivity (p) as a function of temper- 
ature in the heavy fermion compound CeCug. 




FIG. 2. The NMR relaxation rate (1/Ti) and electrical resistivity in the ab-plane (pab) ^ a 
function of temperature in the normal state of the high Tc compound Lai.85Sro.i5Cu04. 
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FIG. 3. Impurity configurations and energy levels in (a) the symmetric Anderson model and 
(b) the asymmetric Anderson model. 
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FIG. 4. The generalized Anderson model as a three-level system. The wavy lines represent the 
conduction electron bath. 
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FIG. 5. A typical hopping sequence in the atomic representation along the imaginary time 
axis T G [0,(3 = l/T]. Ti, for i = l,...,n, labels the time at which the impurity hops from one 
configuration to another. 
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FIG. 6. Schematic picture showing that the fugacities of the Coulomb gas representation cor- 
respond to the dimensionless quantum transition amphtudes between the impurity configurations. 
Ut = tCo is the charge fugacity, and yj = J±S,o the spin fugacity. is the inverse energy cutofi'. 
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FIG. 7. Phase diagram of the generahzed Anderson model Eq. (|l|) in the mixed valence regime. 
Here et and ej label the charge and spin stiffness constants defined in the text. The vertical thick 
line, the horizontal thick line, and the dashed line are the boundaries between the different mixed 
valence states. The dashed line is schematic. 
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FIG. 8. Hopping sequences in the atomic representation for (a) the usual spin-Kondo problem 
and (b) the charge-Kondo problem, i.e., the resonant- level model. 
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FIG. 9. The RG flows for (a) the usual spin-Kondo problem and (b) the resonant-level model. 
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FIG. 10. Crossover diagrams in terms of temperature (T) and the impurity-level (E^) (a) for the 
strong coupling phase where A* is the renormalized resonance width; and (b) for the intermediate 
phase where labels the critical impurity level. 
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FIG. 11. Schematic picture of the lattice model and its reduction to an effective impurity model 
with self-consistent conduction electron bath in the large D limit. 
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FIG. 12. Numerical results of the d— and c— electron Green's functions vs. the Matsubara 
frequency ujn for a set of parameters for which the solution is a non-Fermi liquid metallic state. 
Solid lines come from the self-consistent exact diagonalization method and the solid squares are 
from quantum Monte Carlo with /3 = 64. 
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FIG. 13. (a)Crossover in terms of temperature (T) and chemical potential (/i) for the inter- 
mediate phase; (b) Electron density (n) versus chemical potential for the intermediate phase. 
~ where > 1 is the renormalized charge stiffness constant. 




critical point. The dashed line is schematic. The phases are described in the text. 
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